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Abstract 

We derive the form of the Belavkin-Kushner-Stratonovich equation describing the filtering of a continuous ob- 
served quantum system via non-demolition measurements when the statistics of the input field used for the indirect 
measurement are in a general coherent state. 

Dedicated to Robin Hudson on the occasion his 70th birthday. 



cn 1 Introduction 

' (— I 'One of the most remarkable consequences of Hudson-Parthasarathy quantum stochastic calculus [3D] is Belavkin's for- 
O limitation of a quantum theory of filtering based on non-demolition measurements of an output field that has interacted 
with a given system [H [5j [6j [7]. Specifically, we must measure a particular feature of the field, for instance a field 
quadrature or the count of the field quanta, and this determines a self-commuting, therefore essentially classical, stochas- 
tic process. The resulting equations have structural similarities with the classical analogues appearing in the work of 
L_jKallianpur, Striebel, Kusnher, Stratonovich, Zakai, Duncan and Mortensen on nonlinear filtering, see [HI [551 1231 [57| ■ AS 
ji consequence, earlier models of repeated quantum photon counting measurements developed by Davies [131 114j could 
^ be realized through a concrete theory: this was first shown by taking the pure-jump process limit of diffusive quantum 
qq filtering problems [3|. 

t> ; 

There has been recent interest amongst the physics community in quantum filtering as an applied technique in quantum 
• feedback and control (TJ [2j [9l [TOJ [11] [TBI fUl EJ [25l [26] . An additional driver is the desire to go beyond the situation 
6f a vacuum field and derive the filter for other physically important states such as thermal, squeezed, single photon 
states, etc. In this note we wish to present the filter for non-demolition quadrature and photon-counting measurements 
i~ I jvhen the choice of state for the input field is a coherent state with intensity function f3. The resulting filters are a 
deformation of the vacuum filters and reduce to the latter when we take j3 = 0, this is perhaps to be expected given 
that the coherent states have a continuous-in-time tensor product factorization property. We derive the filters using the 
reference probability approach, as well as the characteristic function approach. 
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1.1 Classical Non-linear Filtering 

We consider a state based model where the state X t evolves according to a stochastic dynamics and we make noisy 
observations Y t on the state. The dynamics-observations equations are the SDEs 

dX t = v(X t )dt + a x (X t )dWr C , (1) 
dY t = h{X t )dt + a Y dW° hs (2) 

and we assume that the process noise W ploc and the observation noise W b s are uncorrelated multi-dimensional Wiener 
processes. The generator of the state diffusion is then 



1 



where Sxjf = (Tx&x- The aim of filtering theory to obtain a least squares estimate for the state dynamics. More 
specifically, for any suitable function / of the state, we would like to evaluate the conditional expectation 



7r t (/) :=E[/(X t )|7f], 

t] 

1.1.1 Kallianpur-Striebel Formula 



with IFK being the er-algcbra generated by the observations up to time t. 



By introducing the Kallianpur-Striebel likelihood function 



L t (x|y) = exp J jft, 



(x s ) T dy s - -h (x s ) T h (x s ) ds 



for sample state path x = {x s : < s < t} conditional on a given sample observation y = {y s : < s < t}, we may 
represent the conditional expectation as 



f csm f(xt)L t (*\y)V[d> 
[J) = - 



fc S [o,t] L * Wy) F [ dx ] 

where P is canonical Wiener measure and 



°t (/) 
o~t (1) 



^t(/)M= / /(x t )L*(x|FH)P[dx]. 

JC*[0,t] 

1.1.2 Duncan-Mortensen-Zakai and Kushner-Stratonovich Equations 

Using the Ito calculus, we may obtain the Duncan-Mortensen-Zakai equation for the un-normalized filter a t (/), and the 
Kushner-Stratonovich equation for the normalized version w t (/). These are 

da t (/) = a t (£/) dt + a t (fh T ) dY t , 

dirt (/) = TT t (£/) dt + [iT t {.fh T ) - 7r t (/) 7T t (/i T )] df t) 

where (It) are the innovations: 

d/ t := dY t - TTt (/i) dt, I (0) = 0. 

We note that there exist variants of these equations for more general processes than diffusions (in particular for point 
processes which will be of relevance for photon counting) , and for the case where the process and observation noises are 
correlated. 



1.1.3 Pure versus Hybrid Filtering Problems 

We remark that we follow the traditional approach of adding direct Wiener noise W bs to the observations. We could of 
course consider a more general relation of the form dY t = h (X t ) dt + aydW° hs but for constant coefficients ay a simple 
rearrangement returns us to the above setup. 

The situation where we envisage dY t = h (X t ) dt + ay (X t ) dW° hs , with ay a known function of the unobserved state, 
must be considered as being too good to be true since we can then obtain information about the unobserved state by just 
examining the quadratic variation of the observations process, since we then have dY (t) dY (t) = ay(X t )ay (X t ) dt. 
For instance, in the case of scalar processes, if we have ay (X) = "y\X\ then knowledge of the quadratic variation yields 
the magnitude \X t \ of the signal without any need for filtering. Such situations rarely if ever arise in practice, and one 
naturally restricts to pure filtering problems. 



2 Quantum Filtering 



Wc wish to describe the quantum mechanical analogue of the classical filtering problem. To begin with, we note that 
in quantum theory the physical degrees of freedom are modeled as observables, that is self-adjoint operators on a fixed 
Hilbert space t). The observables will generally not commute with each other. In place of the classical notion of a state, 
we will have a normalized vector ip G f) and the averaged of an observable X will be give by the real number (ip\Xip). 
(Here we following the physicist convention of taking the inner product to be linear in the second argument </> and 

conjugate linear in the first ip.) More generally we define a quantum state to be a positive, normalized linear functional 
E on the set of operators. Every such expectation may be written as 

E [X] = tr h {gX} 

where g is a positive trace-class operator normalized so that tv^g = 1. The operator g is referred to as a density matrix. 
The set of all states is a convex set whose extreme points correspond to the density matrices that are rank-one projectors 
onto the subspace spanned by a unit vectors if> G f). 

To make a full analogy with classical theory, we should exploit the mathematical framework of quantum probability 
which gives the appropriate generalization of probability theory and stochastic processes to the quantum setting. The 
standard setting is in terms of a von Neumann algebra of observables over a fixed Hilbert space, which will generalize the 
notion of an algebra of bounded random variables, and take the state to be an expectation functional which is continuous 
in the normal topology. The latter condition is equivalent to the cr-finiteness assumption in probability theory and results 
in all the states of interest being equivalent to a density matrix. 

In any given experiment, we may only measure commuting observables. Quantum estimation theory requires that the 
only observables that we may estimate based on a particular experiment are those which commute with the measured 
observables. In practice, we do not measure a quantum system directly, but apply an input field and measure a component 
of the output field. The input field results in a open dynamics for the system while measurement of the output ensures 
that we met so called non-demolition conditions which guarantee that quantum measurement process itself docs not 
destroy the statistical features which we would like to infer. We will now describe these elements in more detail below. 

2.1 Quantum Estimation 

We shall now describe the reference probability approach to quantum filtering. Most of our conventions following the 
presentation of Bouten and van Handel [TT] . 

Let 21 be a von Neumann algebra and E be a normal state. In a given experiment one may only measure a set of 
commuting observables {Y a : a £ A}. Define the measurement algebra to be the commutative von Neumann algebra 
generated by the chosen observables 

ffl = vN{r Q : a G A} C 21. 
We may estimate an observable X G 21 from an experiment with measurement algebra DJl if and only if 

X G M' := {A G 21 : [A,Y] =0,V7e M}, 

That is, if it is physically possible to measure X in addition to all the Y a . Therefore the algebra vN{X, Y a : a G A} 
must again be commutative. We may then set about defining the conditional expectation of estimable observables onto 
the measurement algebra. 

Definition 1 For commutative von Neumann algebra DJt, the conditional expectation onto DJl is the map 

E[- | 971] : M' H- m 

by 

E[E[X \M}Y}= E[XY], \/Y G 9Tt. (3) 



In contrast to the general situation regarding conditional expectations in the non-commutative setting of von Neumann 
algebras |24| . this particular definition is always nontrivial insofar as existence is guaranteed. Introducing the norm 
|^4|| 2 := E[A^A], wc see that the conditional expectation always exists and is unique up to norm-zero terms. It moreover 
satisfies the least squares property 

\\X-E[X | 371] || < \\X -Y\\, VFeOn. 

As the set vN{X, Y a : a e A} is a commutative von Neumann algebra for each X £ 971', it will be isomorphic to the 
space of bounded functions on a measurable space by Gelfand's theorem. The state induces a probability measure on 
this space and we may obtain the standard conditional expectation of the random variable corresponding to X onto the 
er-algebra generated by the functions corresponding to the Y a . This classical conditional expectation then corresponds 
to a unique element E[X | 971] £ 971 and this gives the construction of the quantum conditional expectation. 

We sketch the conditional expectation in figure 1. Note that while this may seem trivial at first sight, it should be stressed 
that the commutant 971' itself will typically be a non-commutative algebra, so that while our measured observablcs 
commute, and what we wish to estimate must commute with our measured observables, the object we can estimate need 
not commute amongst themselves. 

/ — 7, — r 
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Figure 1: Quantum Conditional Expectation E[-|97t] 
The following two lemmas will be used extensively, see [TT] and (T2"] . 

Lemma 2 (Unitary rotations) Let U be unitary and define E [X] := E [t/^X/j] and let 971 = WVMJ . Then 

E[U f XU | 971] = JJ+E [X | 971] U. 

Here we think of going from the Schrodinger picture where the state E is fixed and observables evolve to U*XU, to the 
Heisenberg picture where the state evolves to E and the observables are fixed. Lemma [2] tells us how we may transform 
the conditional expectation between these two pictures. 

Lemma 3 (Quantum Bayes' formula) Let F <E 

E F [X 

Proof. For all Y G 971, 

E [E [F^XF | 971] Y] = E [F^XFY] 

= E F {XY], since [F, Y] = 0, 

= E F [E F [X | 371] y] 

= E [F^FEp [X | 371] Y] , since F e 971' 

= E [E [F^F | 971] E F [X \Tt)Y]. 



( 



371' with E [F^F] = I 
[F^XF | 37t] 



id set . 



[X] := E [F^XF] . Then 



371] 



E [FiF | 37t] 



Note that the proof only works if F £ 371'! 



2.2 Quantum Stochastic Processes 

Wc begin by reviewing the theory of quantum stochastic calculus developed by Hudson and Parthasarathy [20] which 
gives the mathematical framework with which to generalize the notions of the classical Ito integration theory. 

We take M. + = [0, oo). Wc shall denote by -L 2 ymm (R™ ) the space of all square-integrable functions of n positive variables 
that are completely symmetric: that is, invariant under interchange of any pair of its arguments. The Bose Fock space 
over L 2 (M. + ) is then the infinite direct sum Hilbcrt space 

oo 

■8 := ^symmP^™) 
n=0 

with the n = space identified with C. An element of $ is then a sequence W = (VvO^Lo with ip n £ m (M") and 
ll^ll = J2^Lo Jjo oo)" l^n ' ' ' ' I 2 dti • • • dt n < oo. Moreover, the Fock space has inner product 



oo - 

* I $ ) =Y. / tfn fa,"' > *»)>(*!. ■" ,tn). 

ni[0,Oo)» 



Th physical interpretation is that = (VvO^Lo ^ $ describes the state of a quantum field consisting of an indefinite 
number of indistinguishable (Boson) particles on the half-line R + . A simple example is the vacuum vector defined by 

n-.= (i,o, o,-") 

clearly corresponding to no particles. (Note that the no-particle state is a genuine physical state of the field and is not 
just the zero vector of J!) An important class of vectors are the coherent states \1/ defined by 

[*C8)] n (ti,- ) tn):=e-ll"ll !, -i=)9(t 1 )-..i8(t n ), 

vn! 



for (i £ i 2 (R + ). (The n — component understood as e I'* 3 " .) The vacuum then corresponds to VP (0). 
For each t > we define the operators of annihilation B (t) , creation B* (t) and gauge A (t) by 



Jo 

1 - 

[B* (*!,■■• ,*») := -^X)lp,,t](*i)^n-l(tl,- 

3=1 



' 3 ' 



[A(t)*] n (ti,-.. ,t n ) := YlModfaWn {h,--- ,t n ). 

3 = 1 

The creation and annihilation process are adjoint to each other and the gauge is self-adjoint. Wc may define a field 
quadrature by 

Q{t) = B{t)+B* (t) 

and this yields a quantum stochastic process which is essentially classical in the sense that it is self-adjoint and self- 
commuting, that is [Q (t) , Q (s)] = for all t, s > 0. We remark that for each 6 £ [0, 2ix) we may define quadratures 
Qe (t) = e~ %B B (t) + e l6 B* (t) which again yield essentially classical processes, however different quadratures will not 
commute! For the choice of the vacuum state, {Q (t) : t > 0} then yields a representation of the Wiener process: for real 
fc(.) 

We also note that {A (t) : t > 0} is also an essentially classical process and for the choice of a coherent state yields a 
non-homogeneous Poisson process: for real k (■) 



(* 08) | e 1 ^ & « dA W* [fi)) = exp f°° \/3 (t) | 2 (e ww - l) 



di. 



Table 1: Quantum Ito Table 
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We consider a quantum mechanical system with Hilbert space f) being driven by an external quantum field input, the 
quantum field will be modeled as an idealized Bose field with Hilbert space T [L 2 (R+, dt)) which is the Fock space 
over the one- particle space L 2 (M. + ,dt). Elements of the Fock space may be thought of as vectors ^ = (B^ =0 i[>n where 
ipn — ipn {ti, ■ ■ ■ , t n ) is a completely symmetric functions with 



OO r. 

^2 IV>« (hr ■ ■ ,t n ) \ 2 dti ■ ■ -dt n < oo. 



n=Q J [0,00 

The Hudson-Parthasarathy theory of quantum stochastic calculus gives a generalization of the Ito theory of integration 
to construct integral processes with respect to the processes of annihilation, creation, gauge and, of course, time. This 
leads to the quantum Ito table [TJ 

We remark that the Fock space carries a natural filtration in time obtained from the decomposition ^ = $ t ] ® ^ into 
past and future subspaces: these are the Fock spaces over L 2 [0,t] and L 2 (t, oo) respectively. 

2.3 Continuous-Time Quantum Stochastic Evolutions 

On the joint space f) <8> -J, we consider the quantum stochastic process V(-) satisfying the QSDE 



(S-I)® dA(t) +L®dB* it) - L*S ® dB (t) - {\l*L + iH) ® dtX V 



with V(0) = 1, and where S is unitary, L is bounded and H self-adjoint. This specific form of QSDE may be termed 
the Hudson-Parthasarathy equation as the algebraic conditions on the coefficients are necessary and sufficient to ensure 
unitarity (though the restriction for L to be bounded can be lifted). The process is also adapted in the sense that for 
each t > 0, V (t) acts non-trivially on the component f) £g> 3t] and trivially on $u. 

2.3.1 The Heisenberg-Langevin Equations 

For a fixed system operator X we set 

jt (X) := ft {t) [ X ®I]V (t) . (4) 

Then from the quantum Ito calculus we get 

djt (X) = it (£ U X) ® dA (t) + 3t (C W X) ® dB* (t)) + Jt {CoiX) ® dB (t) + j t (£ 00 X) ® dt (5) 

where the Evans-Hudson maps £ M „ are explicitly given by 

CnX = S*XS-X, 
C W X = S*[X,L], 
j£qiX = [L* : X]S 

CooX = £{L,H) 

and in particular Cqq takes the generic form of a Lindblad generator: 

C [l ,b) = h*[X,L] + ±[L*,X}L-t[X 1 H}. (6) 
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2.3.2 Output Processes 

We introduce the processes 



Figure 2: Input - Output component 



B out (t) 
A out (t) 



yt (t) [I (8 B (t)] V (t) , 
V* (i) [J <g> A (t)] V (t) . 



(7) 



We note that we equivalently have i? out (i) = V* (T) [1 <g> i? (t)] y (T), for t <T. Again using the quantum Ito rules, we 
sec that 



dBout = jt(S)dB(t) + j t (L)dt, 

dA out = dA(t)+j t (L*S)dB*(t)+jt(S*L)dB(t)+jt(L*L)dt. 



(8) 



2.3.3 The Measurement Algebra 

We wish to consider the problem of continuously measuring a quantum stochastic process associated with the output 
field, we shall chose to measure an observable process of the form 



Y out (t) := V(t) + [I® Y in (t)}V(t) 
which corresponds to a quadrature of the field when 

Y in (t)=Q(t) = B*(t)+B(t), 
or counting the number of output photons when 

y in (<) = A(t). 

We introduce von Neumann algebra 

2}^vN{y in ( S ):0< S <i}, 
and define the measurement algebra up to time t to be 

2)° ut = vN {y out (s) : < s < t) = V (i) f ifffiV (i) 

Note that both algebras are commutative: 



(9) 



(10) 



for T = t V s. 



[y out (t) , y out ( a )] = v (r) f (/ gi [y in (t) , y in (*)] ) v (t) = o 

The family |?)°] Ut : t > o| then forms an increasing family (filtration) of von Neumann algebras. 



2.3.4 The Non-Demolition Property 

The system observables may be estimated from the current measurement algebra 

j t (X)e (2)?] ut )'. 
The proof follows from the observation that for t > s 

[j t (x) , y° ut ( s )] - y (rf [x ® i, i ® y in ( s )] y (t) = o. 



(ii) 



2.4 Constructing The Quantum Filter 

The filtered estimate for j t (X) given the measurements of the output field is then 

tt 4 (X) := E [j t (X) | 2)°, u 
Let Ef [X] = E [j t (X)], then by lemma [2] 

tt, (x) = e [ Jt (x) | = y (t)+ £ t [x|2j^] y (t) . 

2.4.1 Reference Probability Approach 

Suppose that there is an adapted process F (•) such that F (t) G and E t [X\] = E F (t^ (X <g) 1) F (t) 

system operators X, then by lemma [3] 



for all 



3 



PO I ?)°] u 



V(tyE t [X\%™ V(t) 
E I"f 



E 


F (i) f (X <g> 1) F (t) | 




E 







7(t) 



This is essentially a non-commutative version of the Girsanov transformation from stochastic analysis. The essential 
feature is that the transformation operators F (t) giving the change of representation for the expectation lie in the 
commutant of the measurement algebra up to time t. 



We therefore obtain an operator-valued Kallianpur-Striebel relation 

o-t (X) 



°t (1) 



which may be called the quantum Kallianpur-Striebel, where 

a t (X) := V (tf E \f (tf (X ® 1) F (t) \ 2)^1 V (t) 



(12) 



(13) 



3 Coherent State Filters 

We shall consider the class of states 
of the form 



E' 3 [.] = (^|.^> 

V/ = <j> ® # (/?) 



(14) 
(15) 



where is a normalized vector in the system Hilbert space and VP (/3) is the coherent state with test function j3 G L 2 [Q, oo). 
We note that 

dB (t) V(P) = (t) dt (J3) , 
dA(t) *C8) - P(t)dB*(t)*(fi). 



We see that 



where 



E*[<& (X)]=E^[jt (cMx)]dt 
C m X = CooX + /3 (t)* C W X + (3 (t) CoiX + \(3 (t) \ 2 £ n X 



(16) 
(17) 



The generator is again of Lindblad form and in particular we have 

rW) = C, 

^ — J ~(L?( t '>,HP'. t )) 

with 

If® =S0(t)+L, H ® = j. (L + j9(t) - Lj3(t)*) . (18) 

We may define a parametrized family of density matrices on the system by setting trj, {g t X} = E^[jt (X)], in which case 
we deduce the master equation 

where the adjoint is defined through the duality trj, {g£X} =tr^ {C'gX}. 
From the input-output relation for the field 

cLBout = jt (S) dB (t) + j t (L) dt 

we obtain the average 

E^ [dB out ] = {j t (S) 13 (t) + j t (L)} dt = j t (if®) dt. 



3.1 Quadrature Measurement 

We take Y ln (t) = B (t) + B* (t) which is a quadrature of the input field. Setting tp (t) = V (t) tp we have that 



dip (t) = [(S -l)/3 (t) + L] dB* (t) i/> (t) 



L*SP(t) + -L*L + iH 



dtijj(t) (19) 



At this stage we apply a trick which is essentially a quantum Girsanov transformation. This trick is due to Belavkin [5] 
and Holevo [TH] ■ We now add a term proportional to dB (t) ip (t) to get 

dip (t) = ((5 -1)13 (t) + L) [dB* (t) + dB (t)} ip (t) 

L*Sf3 (t) + ^L* L + iH + ((S - 1) (3 (t) + L) (3 (i) j dtip(t) 



L t dY' ln {t)ip(t) + K t dt^{t). 



where 



L t = L + (S-I)(3(t) = L^ -f3(t), 

K t = -L*S(3{t)--L*L-iH-L f) ®l3{t)+l3{tf . 

It follows that ip (t) = F (t) tp where F (t) is the adapted process satisfying the QSDE 

dF (t) = L t dY in (t) F (i) + K t dt F (t) , F (0) = I. 
Moreover F (t) is in the commutant of 2}™ and therefore allows us to perform the non-commutative Girsanov trick. 
From the quantum Ito product rule we then see that 

d [F* (t) XF (t)} = F* (t) (xlt + L* t X^j F (t) dY in (t) + F* (t) (L* t XL t + XK t + K t X} F (t) dt 
and this leads to the SDE for the un-normalizcd filter 

da t (X) = a t (xL t + L* t X^j dY out (t) + a t (l*XI t + XK t + K t X \ dt. 
After a small bit of algebra, this may be written in the form 

da t {X) = a t (xL t + L* t X^j [dY out (t) - ((3 (t) + /3 (*)*) dt] + a t (c m X^ dt. 



This is the quantum Zakai equation for the filter based on continuous measurement of the output field quadrature. 



To obtain the quantum Kushncr-Stratonovich equation we first observe that the normalization satisfies the SDE 

da t (I) = a t {i t + L*) [dY mt (t) - (/3 (t) + (3 (t)*) dt] 

and by Ito's formula 

2 



i a t (L t + L*] at(L t + L*)' 

i— = V ' [dY^ (t) - (ft (t) + p dt] + V > 

t ) a t (I) a t (I) 



The product rule then allows us to determine the SDE for irt (X) = 

°t (I) 

dirt (X) = n t (c m X^ dt + [ir t (xL t + L* t x) - n t (X) n t (l t + L*) } dl (t) , 
where the innovations process satisfies 



dl (t) = dY out (t) - [ n (L t + Lfj +/3(t)+j3 (t)* 



dt. 



We note that the innovations is martingale with respect to filtration generated by the output process for the choice of 
probability measure determined by the coherent state. 

3.1.1 The Quadrature Measurement Filter for a Coherent state 

In it convenient to write the filtering equations in terms of the operators L^W . The result is the Belavkin-Kushner- 
Stratonovich equation for the filtered estimate based on optimal estimation of continuous non-demolition field-quadrature 
measurements in a coherent state /?(•). 



d-K t 

with the innovations 



(X) = 7T t (c m x) dt + {n t (xL™ + L^>X\ - ir t (X) ir t (if® + L^>) } dl^ ad (t) , 

d/qU ad ^ = dY out ^ _ nt r L m + L m*^ dL (20) 



3.2 Photon Counting Measurement 

For convenience we shall derive the filter based on measuring the number of output photons under the assumption that the 
function j3 is bounded away from zero. We discuss how this restriction can be removed later. We now set Y ln (t) — A(t). 
We again seek to construct an adapted process F (t) in the commutant of 2)J? such that ip (t) = F (t) ip. We start with 
[T9l again but now note that 

dB* (t)4(t) = -LdA(t)1>{t) 

and making this substitution gives 

dip (t) = -l-.L t dY in (t)^ (t) - (^L*L + iH + L*Sf3 {t)j dt^(t) . 
We are then lead to the Zakai equation 

da t (X) = j^yp^t (l;XL t + m*XL t + L* t Xf3 (t)) dY mt (t) 

~ a t (t^XL*L + ^L*LX + i [X, H] - XL*Sf3 (t) - (3 (t)* S*LX^j dt 

which may be rearranged as 

da t (X) = a t (C^X) dt + TgT^p^ (l* t XL t + P{t)*XL t + L* t Xp (i)) [dY out (t) - \0 (t) \ 2 dt] . 



To determine the normalized filter, we note that 

(i) = T^vi°t (L;L t + p(tyL t + l*j (t)) [dr out (t) - \p (t) \ 2 dt] 



and that 



1/3 (*) 

a t (LtLt + p(t)*Lt + LtP(t) 



dt 



°t (I) a t (i f 

a t (L* t L t +j3{t)*L t + L* t P{t) 



dY out 0) 



a t (I) [a t (LfLt + p(t)*Lt + l\$ (i)) + \0 (t) | s 

An application of the Ito product formula then yields the quantum analogue of the Kushner-Stratonovich equation for 
the normalized filter; 

d-Kt (X) = 7T t (C X) dt + r dl (t) 

n t (L* t L t + f3(t)*Lt + L* t (3 (*)J + \(3 (t) |a 

x {n t (L* t XL t + (3(t)*XL t + L* t Xp (t)) - n (X) n t (L* t L t + p(t)*L t + L*/3 (*))} 

and the innovations process is now 



dl (t) - dY ont (t) - [n t (i* t L t + P(t)*L t + L*j3 (t)) + \0 (t) f 



dt. 



We note that the innovations is again a martingale with respect to filtration generated by the output process for the 
choice of probability measure determined by the coherent state. 

The derivation above relied on the assumption that (3 (t) ^ 0, however this is not actually essential. In the case of a 
vacuum input, it is possible to apply an additional rotation W (t) satisfying dW (t) = [z*dB (t) — zdB* (t) — ^\z\ 2 dt]W (t) 
, with W (0) = /, and apply the reference probability technique to the von Neumann algebra generated by N (t) = 
W (t) A (t) W (t)*: this leads to a Zakai equation that explicitly depends on the choice of z G C however the Kushner- 
Stratonovich equation for the normalized filter will be ^-independent. Similarly for the general coherent state considered 
here, we could take z to be a function oft in which case we must chose z so that f3 (t)+z (t) ^ 0. The Kushner-Stratonovich 
equation obtained will then be identical to what we have just derived. 

3.2.1 Photon Counting Measurement in a Coherent State 

Again it is convenient to re-express the filter in term of L^^' . We now obtain the Belavkin-Kushncr-Stratonovich equation 
for the filtered estimate based on optimal estimation of continuous non-demolition field-quanta number measurements in 
a coherent state /?(•). 

dnt (x) = n (cwx) dt + 1 ^Z>Lm) - "* dinum ^ ' < 21 > 

with the innovations 

dl num (t) = dY out (i) - TT t (L^>L^ dt. (22) 



4 Characteristic Function Approach 

As an alternative to the reference probability approach, we apply a method based on introducing a process C (t) satisfying 
the QSDE 

dC(t) = f(t)C(t)dY(t), (23) 

with initial condition C (0) = /. Here we assume that / is integrable, but otherwise arbitrary. The technique is to make 
an ansatz of the form 

dnt (X) = F t (X) dt + H t (X) dY (t) (24) 



where we assume that the processes Tt (X) and Ht (X) are adapted and lie in 2) t i. These coefficients may be deduced 
from the identity 

E[( n (x)-j t (X))C(t)] = o 

which is valid since C (i) <E 2)t] . We note that the Ito product rule implies I + II + III = where 



/ = E [(d7T t (X) - djt (X)) C (t)} , 
II = E[(Tr t (X)-j t {X))dC{t)], 
III = E [(dirt (X) - dj t (X)) dC (t)} . 

We illustrate how this works in the case of quadrature and photon counting in a coherent state. For convenience of 
notation we shall write St for j t (S), etc. 

4.1 Quadrature Measurement 

Here we have 

dY (t) = S t dB (t) + S^dB (t)* + (L t + L* t ) dt 



so that 



/ = Ep[T t {X)C{t) + Ht(X)(Stl3t + SZl3t + L t +L* t )C{t)}dt 



-E 



dt, 



p [{(c 00 x) t + (c m x) t p t + (c w x) t p* + (£ u x\ \p t \ 2 } c (t) 

II = Ep[{Tr t (X)-X t )f{t)C(t)(StPt + S;p:+L t + L* t ))dt, 
III = Ep[{U t (X)-(C 01 X) t Stfit-(C 11 X) t StPt}f{t)C{t)]dt. 

Now from the identity I + II + III = we may extract separately the coefficients of / (t) C it) and C it) as / (t) was 
arbitrary to deduce 

n t ((n t (X) - X t ) (S t Pt + S* t P* t +L t + L*)) + n t (H t (X) - (C 01 X) t S* t P* t - (£ n X) t S?ft) = 0, 

*t(rt(x) + Ht(x)(StPt + s;p; + L t + L;)- = o. 

Using the projective property of the conditional expectation (7r t o 7r t = 7r t ) and the assumption that Tt (X) and Jit (X) 
lie in 2) t ], we find after a little algebra that 

H t {X) = n t ^XL^+L^*X^-Tr t (X)n t ^L^ +1^*^ 

T t {X) = n(c m x)-n t (X)7r t (L^+L^*), 

so that the equation reads as 

dirt (X) = nt (c m X^j dt + H t (X) [dY (t) - 7T t (L m + ^ (t) *) dt 

4.2 Photon Counting Measurement 

We now have 



dY (t) = dA (t) + L* t S t dB (t) + S* t L t dB (t)* + L*L t dt 



so that 



I 
II 



Ep 
Ep 



{ji (X)+H t (X) (|/3 t | 2 + L* t S t f3 t + S* t L t Pt+L* t L t )}c{t)\ dtEp [{- } C (t) 

(n t (X) - X t ) f (t) C (t) (|/3 t | 2 + L*S t p t + S* t L t p* t + L* t L t 



dt 



III = E 
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H t (X) [\Pt\* + L\S t Pt + StUPl + L\Lt) f (t) C (t) 
l3 [{(CoiX) t S* t L t + (CaiX) t Pt] f (t) C (i)] dt 



dt 



-E, 



{(£uX) t \p t \ 2 + (£ n X) t S* t L t pl}f (t) C (t) 



dt. 



This time, the identity I + II + III = implies 

n t ((n t (X) X t ) L™*L™ + H t (X) if «*if - n t ((C ai X) t S* t L^ + (C u X) t S* t L^fi) = 0, 



n t ^ t (x) + n t (x)Lf ) *L^ t) -[c^x) t ) = 0. 



Again, after a little algebra, we find that 



T t {X) = ■K t (L^x)-Ut{X)-K t (L^>L^), 
so that the equation (J2U) reads as 

dm (X) = n t (CW>X\ dt + Ht (X) \dY (t) - n t (l^>L^ 



dt 



In both cases, the form of the filter is identical to what we found using the reference probability approach. 

5 Conclusion 

Both the quadrature filter ([20]) and the photon counting filter (|2ip take on the same form as in to the vacuum case and 
of course reduce to these filters when we set /3 = 0. In both cases it is clear that averaging over the output gives 

E^ 3 [dnt (X)] = E* 3 \j t X 



dt 

which is clearly the correct unconditioned dynamics in agreement with (|16|) . and we obtain the correct master equation. 

It is worth commenting on the fact that the pair of equations now replacing the dynamical and observation relations 
(|1I2[) arc the Heiscnbcrg-Langcvin equation (0 and the appropriate component of the input-output relation ((8]). The 
process and observation noise have the same origin however the nature of the quantum filtering based on a non-demolition 
measurement scheme results in a set of equations that resemble the uncorrelated classical Kuslmer-Stratonovich equations. 

5.1 Is Quantum Filtering still a Pure Filtering Problem? 

The form of the input-output relations ([8J might suggest that it is possible to learn something about the system dynamics 
by examining the quadratic variation of the output process, however, this is not the case! We in fact have an enforced 
"too good to be true" situation here as the output fields satisfy the same canonical commutation relations as the inputs 
with the result that the quantum Ito table for the output processes -B ut; B* nt and A out has precisely the same structure 
as table [T] Therefore we always deal with a pure filtering theory in the quantum models considered here. 
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